Light Qubit Storage and Retrieval using Macroscopic Atomic Ensembles 
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We present an experimentally feasible protocol for the complete storage and retrieval of arbitrary 
light states in an atomic quantum memory using the well-established Faraday interaction between 
light and matter. Our protocol relies on multiple passages of a single light pulse through the atomic 
ensemble without the impractical requirement of kilometer long delay lines between the passages. 
Furthermore, we introduce a time dependent interaction strength which enables storage and retrieval 
of states with arbitrary pulse shapes. The fidelity approaches unity exponentially without squeezed 
or entangled initial states, as illustrated by explicit calculations for a photonic qubit. 



Faithful storage and retrieval of an unknown quan- 
tum state of traveling light pulses plays an important 
role in quantum information protocols, such as long dis- 
tance communication 0, Q and quantum computation 
■ Several theoretical proposals for deterministic quan- 
tum memory for light have been put forward using sin- 
gle atoms coupled to high finesse cavities j4|, as well as 
atomic ensembles in free space coupled to light 0, IE • 
Recently 0, the first quantum storage of weak coherent 
light pulses with a higher fidelity than achievable clas- 
sically was demonstrated using the off-resonant Faraday 
rotation of light passing through a macroscopic atomic 
ensemble. The subsequent experimental challenge is to 
demonstrate retrieval of the quantum state back to a 
light pulse following for example the proposal 6j. This 
protocol shows better performance, i.e., higher fidelities, 
than achievable classically, but in order to become truly 
applicable for future quantum information protocols, fi- 
delities close to unity must be achieved. This is possi- 
ble with the protocol, but only if one is able to prepare 
highly squeezed initial states, i.e., spin squeezed atomic 
states and squeezed input light beams for mapping and 
retrieval respectively. The difficulty in producing such 
squeezed states limits the realistically achievable fidelity 
of the quantum memory. In this Letter we present a 
new protocol which can reach unity fidelity without any 
squeezing. The requirement for this protocol to reach 
unit fidelity is sufficient interaction strength, which can 
be met by a sufficient number of atoms and photons tak- 
ing part in the interaction. Unlike previous proposals our 
theoretical analysis considers a long light pulse which is 
reflected and redirected so that it is on its way through 
the sample along different directions at the same time, 
see Fig. Q] This dramatically changes the dynamics of 
the whole system, and the evolution is no longer gov- 
erned by the Quantum Non Demolition (QND) coupling 
used in the direct mapping protocol Q and considered in 
various atom-light interface protocols discussed in Refs. 
0, 0, . These protocols have assumed sequential pas- 
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sage of the light through the sample which for pulses of 
millisecond duration as required in Q means that de- 
lay lines of hundreds of kilometers would be necessary. 
In addition to presenting a protocol which works for 
long pulses, we demonstrate that it is possible to tailor 
the interaction strength such that the quantum memory 
may store and retrieve incoming light pulse with arbi- 
trary pulse shapes. Finally, we note that, although being 
an effectively continuous variable system a macroscopic 
atomic ensemble can be used to store a single qubit and 
we shall show explicitly that our protocol enables high 
fidelity storage and retrieval of qubit states in the form 
of superposition states of and 1 photons. 

We investigate a setup with a macroscopic atomic spin 
angular momentum oriented along the x-direction |Tlj or, 
equivalently, two samples with opposite macroscopic spin 
placed in a constant magnetic field See Fig. ^ As 
described in [H, the total angular momentum along 
the x-direction, J Xl is effectively a classical parameter 
and the two remaining components can be rescaled to 
form operators, xa = Jy/^fJx and pa = Jz/V~Jx with 
canonical commutation relation [cc^,^] = i- The light 
state is composed of a strong coherent component lin- 
early polarized in the x-direction and the weak quantum 
mechanical signal to be stored in the y-polarization. For 
this system the Stokes vector component S x is a macro- 
scopic classical quantity and we can define the operators 
x l = S y /\/~S x ~ and pl = S z /^fS^ with commutation re- 
lation [& = i. 

As shown in Fig. Q] our protocol begins with a light pas- 
sage along the z-direction. The off-resonant interaction 
Hamiltonian for this process is given by H oc paPl 0- 
After a 90 deg. rotation in xl~Pl space the light is sent 
through the atomic sample along the y-direction, creating 
a H oc xaxl interaction. To describe the effect of this in- 
teraction we first review the case where the two passes of 
the beam occur one after the other ||| . The output oper- 
ators are then given by p™ 1 = (1 — K 2 )p™ — kx™ , x™' = 

xf + K P % , P ° A ut = {l-K 2 ) P \ l -Kxf , X / 1 = X^+Kpf, 

where k is the integrated interaction strength, k 2 is pro- 
portional to the total number of atoms and the total num- 
ber of photons in the pulse [T3. . The protocol both 
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Figure 1: a) Schematic of the multi-pass protocol. The light 
is sent through the atomic sample along the z- and then along 
the y-direction. After the second passage, we may reflect the 
light back onto the sample to improve the performance by 
two further passages of the light beam, b) Schematic of an 
implementation if the atomic system consists of two atomic 
samples with oppositely oriented mean spin in a homogeneous 
magnetic field as in ||. 



maps the light properties on the atoms and the atomic 
properties on the outgoing light pulse. For k = 1, mean 
values are faithfully stored whereas uncanceled contribu- 
tions from x™ or x™ will limit the fidelity of storage or 



retrieval of e.g. a coherent state to 88% |6|]. As discussed 
in the introduction this deficiency can be remedied by 
initially squeezing these variables (e.g. xa — ► xa/Vc, 
Pa — > PA\/e), in which case the fidelity may be increased 
towards 100%. This however, represents an unachiev- 
able limit due to the insurmountable difficulty in creating 
arbitrarily squeezed states. Furthermore, as mentioned 
above, the requirement of hundreds of kilometer long op- 
tical pulses cannot be avoided for the experimental real- 
ization of Ref . . 

To circumvent this, we propose simultaneous passage 
of the light pulses through the atomic medium in two 
perpendicular directions. We treat this case by split- 
ting the incoming light pulse into segments of duration 
r as discussed in [lfl [l3| • The interaction can then be 
treated sequentially using the coarse grained Hamiltonian 
H Tj i — K T PAPL,i and H' Ti = k t xaxla- The crossing of 
beams in the atomic sample gives rise to standing wave 
interference patterns. Since these will either be station- 
ary or move with a speed of c/y/2 in the 45 degree direc- 
tions with respect to the beams, atoms with a thermal 
velocity distribution will effectively average over these in- 
terferences and only experience the traveling wave com- 
ponents. By taking the r — ► limit, differential equa- 
tions are obtained, which can be solved under the physi- 
cally reasonable assumption that the atomic state evolves 
slowly compared to the time it takes for light to travel 
between the two passages [l-ij: 



x A (t)=x A (0) + / dt'kPi n (t') 
Jo 

ft 

p A {t)=p A {0)e-ti> duk2 - / dt'ke-^' dufi2 X^{t') 

Jo 



xr t (t) = xr(t) + k(t) PA (t) 



(1) 



where for future use we have defined a time dependent 
interaction rate k(t) which is characterized by K 2 ot — 
Jo°° K 2 (t)dt (time dependence of h is implicit in the inte- 
grals). Xl(£) and Pi,(t) are the position and momentum 
operators for the light segment arriving at time t with 
commutation relation \Xi,(t), Pt(t')] = i6(t — t'). As 
can be seen the solution turns out to be very similar to 
the simple QND case reviewed earlier. In particular we 
recognize an asymmetry in which only the initial quan- 
tum fluctuations of one of the atomic quadratures are 
damped, which is exactly the reason why squeezing of 
the initial states is required. 

The result in itself is noteworthy, because it gives the 
first indication that retrieval of stored states is experi- 
mentally feasible for the system of Ref. 8j. Indeed we 
have devised a protocol, which uses this two pass retrieval 
in combination with the single pass storage demonstrated 
in [jj to achieve an ideal quantum memory if infinitely 
squeezed states are available |14| ] . Our main purpose here 
is to show that one may devise a protocol which does not 
require squeezing, and we shall not go into more detail 
with this direct mapping 8] /two-pass retrieval ("1+2") 
protocol here, except for a comparison of its performance 
with the main novel protocol of this paper in which the 
light is reflected back through two further passages of 
the atomic medium, see Fig.^ In the resulting four pass 
protocol the light will experience an interaction sequence 
pp, xx, xx, and pp, which effectively restores the sym- 
metry between the atomic quadratures. As in the two 
pass case a differential equation for the time evolution 
of the atomic and light variables can be formulated and 
solved to obtain [l4|: 

i-t 

x A {t) = x A (0)e- 2 f° duk2 + 2 / dt'ke- 2 ti d ^p^(t') 



Pa 



xr\t) = x?{t) +2~K(t) PA (t) 
pr t [t)=p?{t)-2k{t)x A [t). 

By introducing the annihilation operators a = (x + 
ip) I \[2 and A = (X + iP) / \[2 the entire interaction can 
be written as a beam splitter relation: 

pt 

a A (t) = d A (0)e- 2 I° duh2 -2i dt'ke- 2 ^ duk2 A{t'). 

Jo 

(3) 

To quantify the quantum memory performance, we cal- 
culate the combined storage and retrieval fidelity F = 
(*ideai|p|*idcai), where ^ideai is the ideal state and p 
is the retrieved density matrix. We consider the stor- 
age and retrieval of an unknown qubit stored in the zero 
and one photon subspace of the light field, and calculate 
the fidelity averaged over all orientations on the Bloch 
sphere. For the "4+4" protocol, this can be calculated 
directly from Eq. by noting that imperfect opera- 
tion results in admixture of vacuum noise into the mode 



{t) = p A (0)e- 2 -fo duk2 -2 dt'Re- 2 ti d «« 2 X™(i') 
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Figure 2: Combined storage and retrieval fidelity of a qubit 
with constant interaction strength for one pass storage and 
two pass retrieval, e = 4, (crosses) and four pass storage and 
four pass retrieval (full curve). In the table complete storage 
and retrieval fidelities (obtained from the optimum of curves 
as the ones shown) for two different retrieval protocols with 
varying degrees of squeezing are presented along with the op- 
timum values of the accumulated interaction strength, «tot- 



we are interested in, and the effect (or an actual source 
of loss or decoherence) is therefore fairly easily deter- 
mined, e.g., by noting that within the and 1 photon 
subspace (0 corresponds to the relation for a decaying 
two-level system. The calculations for protocols involving 
either a single or a two-pass component are significantly 
more complicated involving orthogonal mode decompo- 
sition and appropriate differentiations of the Husimi Q- 
function and we shall not go into details here 0|. In 
Fig. [21 we show two examples of the behavior of the fi- 
delity as a function of the total interaction strength n 2 ot 
for constant k(t). Experimentally k is proportional to 
the amplitude of the classical light field, and hence the 
total number of photons is given by the time integral of 
k 2 . As discussed below the noise introduced by spon- 
taneous emission each time the light passes through the 
sample is proportional to the number of photons. Hence, 
in order to be able to compare protocols involving dif- 
ferent number of passes we compare their performance 
as a function of total numbers of photon passes which, 
e.g, for four pass storage followed by four pass retrieval is 
proportional to n 2 ot = 8 J dtk 2 (t). The crosses show the 
result of the direct mapping protocol with n 2 ot = 2 and 
retrieval of light with a two pass interaction (" 1+2" ) . 
Both the initial atomic state and the light of the retrieval 
pulse are initially squeezed by a factor e = 4. The solid 
curve shows the fidelity for the " 4+4" protocol. The table 
summarizes the optimal combined storage and retrieval 
fidelities at different squeezing factors for the "1+2" and 
"1+4" protocols. In each case the total accumulated in- 
teraction strength is also given. We note that with a 
moderate amount of squeezing, fidelities well above the 
classical boundary of 2 /3 can be achieved. We also note 
that the "4+4" protocol does not require any squeezing, 
and with a moderate n 2 ot = 4.8 we get F — 88% with 
constant coupling strength. 

As a departure from all previous proposals we propose 
to increase the achievable degree of fidelity and at the 
same time facilitate the storage of light pulses with ar- 
bitrary temporal profile by appropriately varying the in- 



teraction constant throughout the pulse. This is readily 
achieved since the time dependent interaction strength 
k(t) is proportional to the amplitude of the strong x- 
polarized driving field, which can be varied experimen- 
tally using e.g. an electro-optic modulator (EOM). In or- 
der to store an incoming light pulse with a temporal mode 
function f(t) so that it is described by the annihilation 

operator af = ^ dtf(t)Af(t), we want a A (T) = af 
which is equivalent to 

1 ^^-2 1 df 
k dt f dt 

If we are interested in storing the state of a light pulse 
described by any (real) f(t) this differential equation can 
be solved to obtain the appropriate shape of the inter- 
action strength k. Hence by suitable tailoring of k the 
protocol allows for ideal storage (and retrieval) of an in- 
coming light state of arbitrary shape without any initial 
squeezing. Note, that the shape of k for storage and re- 
trieval may be different so that the retrieved pulse may 
be transformed into any desirable shape. As a specific ex- 
ample let us consider a constant field mode of duration 
T with f(t) = 1/VT. In the case of mapping, the initial 
atomic state is damped exponentially but so is the early 
part of the input light pulse. This can be counteracted 
by an increased interaction strength for the front part of 
the pulse, and Eq. J3J gives k(t) — -^-^ for the optimal 
mapping interaction. When retrieving a stored state the 
rear end of the light pulse reads out a damped atomic 
state and the divergence needs to be placed at this end 
of the pulse. 

Another interesting result can be obtained if we let 
k(t) — 2 ^ t+T ) where T is the pulse length. In this case 
we get a 50/50 beam splitter between the light and atoms. 
This adds to the growing toolbox of interesting interac- 
tions between atomic ensembles and light pulses, which 
pave the way for a number of hybrid light-atom protocols, 
e.g., for teleportation and entanglement swapping. 

Our protocol for optimum storage and retrieval of pho- 
tonic states requires divergences in the coupling strength, 
which cannot be achieved experimentally. To quantify 
the effect of finite k, we assume that we truncate k 
by making it constant close to the divergences, 
and 1/y/T — t for the mapping and retrieval respectively, 
e.g., for the retrieval k = min [l/2yJT^t,(j)\ where <j) 
is a chosen constant, see inset of Fig. [3J The fidelity 
for the combined storage and retrieval can be calculated 
analytically 0|, and for large interaction the fidelity ap- 
proaches unity rapidly: 

F _> l _ V£_i£ cxp (_ K 2 ot/2 ) , Ktot » l . ( 5 ) 

In Fig. El we compare the "1+2" and the "4+4" pro- 
tocols, both with time-dependent k. For the "1+2" we 
include curves for which both x™ and xf are squeezed 
with a factor e = 1,4 . We see that fidelities of the or- 
der of 93 % can be achieved with experimentally feasible 



4 




3 4 5 6 7 o B 9 10 11 12 

K ,o, 



Figure 3: The fidelity of mapping and subsequent retrieval of 
a qubit with different truncations of the interaction strength 
@. Curves are shown for the one pass direct mapping fol- 
lowed by two-pass retrieval with and without squeezing, e = 4, 
of the initial states (dashed and dotted respectively) and for 
two truncated four pass pulses (full curve). Inset: example of 
a truncated interaction strength for retrieval. 

squeezing. As expected the four pass protocol without 
squeezing approaches unit fidelity already for experimen- 
tally achievable values of «t ot , whereas significant squeez- 
ing is necessary for the " 1+2" protocol. We thus see that 
the "4+4" protocol solves the challenge of feasible high 
fidelity quantum memory. Note, that this protocol in- 
volves a significant amount of passages so if the atomic 
sample is contained within glass cells as in light re- 
flection loss at each interface would decrease the fidelity. 
This decrease, however, can be reduced by adding anti 
reflection coating to the glass cells or by having the addi- 
tional mirrors within the vacuum chamber in a trapped 
atom setup. 

The effect of spontaneous emission of light by the 
atoms should also be addressed. For a single pass 
through a gas with resonant optical density a the spon- 
taneous emission probability per atom r\ is given by r\ = 
K tot/ a U3- For several passes the probability of sponta- 
neous emission is proportional to the number of passes, 
which is why we multiply the single passage Kjr ot by the 
number of passes when comparing different schemes. Be- 
cause of the spontaneous emission it is important to keep 
Kj ot fairly small. The effect of spontaneous emission will 
reduce the fidelity of the protocol by an amount propor- 
tional to r), so that in combination with Eq. J5J the total 
error (1 — F) is of the order of ~ exp(— arf) + r\. Op- 
timizing this expression with respect to r\ we find that 



the error scales as \og(a)/a. It is thus advantageous to 
make a as large as possible, something which is readily 
achieved in trapped cold samples • Here large optical 
depths are achieved for ultra compact samples of atoms, 
thus preserving the prospect of future miniaturization. 

The present protocols can also be extended to employ 
two atomic memory units (ensembles) A and B to store or 
retrieve a qubit originally represented by a single photon 
in two modes L and M, \iP)lm = Q! |10)la/ + /?|01)lm- 
This is important since this encoding is used in a vast 
majority of current experiments on optical quantum in- 
formation processing with discrete variables. The fidelity 
of storage and/or retrieval can be calculated using similar 
techniques as for the single-mode qubit |l4j . 

In conclusion, existing protocols for using the off- 
resonant Faraday interaction to retrieve an unknown 
light state stored in a macroscopic atomic sample 8] are 
very impractical because they assume sequential passage 
of the light beam which requires prohibitively long delay 
lines. Furthermore, even if the protocols could be real- 
ized experimentally they predict unity fidelity only in the 
limit of infinitely squeezed input states. 

Our solution to both of these problems involves three 
novel features. First of all, we have solved the dynamics 
arising when the probe pulse travels through the atomic 
sample in two orthogonal directions simultaneously. This 
we applied to a novel scheme involving four passages of 
the light. Combined with the third new component: a 
time varying interaction strength, this protocol can be 
used for storage or retrieval symmetrically and works 
with fidelity exponentially approaching unity with in- 
creased interaction strength. No squeezing is required 
and since the interaction strength in the Faraday inter- 
action depends on the number of atoms and the number 
of photons there is no fundamental limit to the achiev- 
able fidelity for this protocol. As a final remark, we note 
that the feasibility of the protocol is accentuated by the 
fact that all basic components are in place already. In 
experiments like 8] one only has to add eight mirrors as 
illustrated in Fig. 
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